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Abstract. We study the lengths of the cycles formed by trajectories in the 
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\ 1. Introduction 

. Bose and Einstein understood 80 years ago that a curious phase transition occurs in a gas 

of non-interacting bosons; it is now commonly refered to as Bose-Einstein condensation. 
Real particles interact, however, and for many years there were doubts that this transition 
takes place in natural systems. London suggested in 1938 that superfluid Helium undergoes 
a Bose-Einstein condensation, and this idea is largely accepted nowadays. Bogolubov 
considered interacting systems; careful approximations allowed him to get back to a non- 
interacting gas, but with a different dispersion relation. See |16j and [S] for more discussion 
Q-r and partial justifications of Bogolubov theory. 

In 1953 Feynman studied the system in the Feynman-Kac representation [3]. The 
partition function can be expanded as a gas of trajectories living in {d + 1) dimensions. 
The extra dimension is commonly refered to as "the time", although it is not related 
to physical time. The situation is illustrated in Fig. ^ A finite system with N particles 
induces a probability on the group SV of permutations of N elements. Feynman considered 



B 

> 

^ ■ the probability for a given particle to belong to a cycle of length n. In the thermodynamic 

limit, there may be strictly positive probability for infinite cycles to be present, and 
Feynman suggested to use this as an order parameter for Bose-Einstein condensation. 

A few years later, in 1956, Penrose and Onsager introduced the concept of "off-diagonal 
long-range order" Formally, it is a correlation between positions x and y given by 
a(x,y) = (c i (x)c(y)}. The system displays off-diagonal long-range order when this corre- 
lation is strictly positive, uniformly in the size of the system and in \x — y\ — > oo. One can 
write a Feynman-Kac version of this correlation, and it involves a special cycle starting at 
x and ending at y; this cycle may wind many times around the imaginary time direction. 
In the limit where x and y are infinitely distant there corresponds a notion of infinite open 
cycle that is reminiscent of Feynman's approach. 

Feynman's order parameter is simpler; it is often used in numerical simulations or 
in order to gain heuristic understanding. On the other hand, everybody agrees that 
Penrose and Onsager order parameter is the correct one. Surprisingly, the question of 
their equivalence is usually eluded, and many physicists implicitely assume equivalence to 
hold. The first mathematical investigation of this question is due to Siito, who showed 
that equivalence holds in the ideal gas. Indeed, he proved that infinite cycles occur in 
the presence of condensation [T2], and that no infinite cycles occur in the absence of 
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Figure 1. The Feynman-Kac representation of the partition function for 
a gas of bosons. The horizontal plane represents the d spatial dimensions, 
and the vertical axis is the imaginary time dimension. The picture shows 
a situation with five particles and two cycles, of respective lengths 4 and 1. 

condensation ^3]; the latter result uses probabilistic methods from the theory of large 
deviations. These results have been extended to mean-field systems in and 

In this paper we explore the links between Feynman cycles and off-diagonal long-range 
order. Let cr(x) denote the off-diagonal correlation between the origin and x E R rf , and 
g(n) denote the density of particles in cycles of length n. We propose the following formula 
that relates both concepts: 

a(x) = ^2 c n {x)g{n) + c 00 (x)g(oo). (1.1) 

n > 1 

Mathematically, the problem is not well posed. Many choices for the coefficients c n are 
possible — a trivial choice is c n (x) = cr(x)/p for all n, including n = oo. We will see, 
however, that there is a natural definition for c n (x) in terms of Wiener trajectories. In 
any case, we conjecture that Eq. (|1.1|) holds with coefficients satisfying 

< c n (x) ^ 1, < Coo (x) 1, 

for all n, x. In addition, we should have 

lim c n (x) = Coo(x) 

71— >OC 

for any x, and 

lim c n (x) = 

|a;|^oo 

for any finite n, but not uniformly in ti\ c.qq{x^) may converge to a strictly positive constant 
c. If c = 1, we get from the dominated convergence theorem that lim^i^oo <r(x) = £>(oo) 
- in which case the off-diagonal long-range order parameter is equal to the density of 
cycles of infinite lengths. 

We establish this formula and these properties in the case of the ideal gas, where we 
show that 

c n (x)= e"* 2 / 4 ^, c OQ (x) = l. (1.2) 

We discuss the validity of the formula (|1.1|) in the interacting gas, proving that these 
properties hold true in a regime without Bose-Einstein condensation. The two order 
parameters should not be always equivalent, however. It is argued in jT3] that they differ 
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when the bosons undergo a regular condensation into a crystalline phase. There is no 
off-diagonal long-range order, but infinite cycles may be present. 

We work in the Feynman-Kac representation of the Bose gas. This representation is 
standard, see e.g. [1] for a clear and concise review, and jS] for a complete introduction. 
We assume the reader to possess some familiarity with it and in Section [21 we directly 
define the main expressions — partition functions, density of cycles, off-diagonal long- 
range order — in terms of space-time trajectories. But basic notions and properties are 
reviewed in Appendix [XJ 

The situation simplifies in absence of interactions; we consider the ideal gas in Section 
|2J where we state and prove the formula that relates the two order parameters. The ideal 
gas is best discussed in the canonical ensemble. Rigorous proofs of macroscopic occupation 
of the zero mode have been proposed and they involve the grand-canonical ensemble, with 
a chemical potential that depends on the volume. Appendix El proposes a simple proof in 
the canonical ensemble. 

Interacting systems constitute a formidable challenge; they are discussed in Section^ 
where partial results are obtained. 

In this paper, we denote finite volume expressions in plain characters, and infinite 
volume expressions in bold characters. Further, we always consider the canonical and 
grand-canonical ensembles where the temperature 1/(3 is fixed; we alleviate the notation 
by omitting the j3 dependence of all quantities. 



2. Feynman cycles and off-diagonal long-range order 

2.1. Partition functions. Our Bose gas occupies a d-dimensional domain D, always a 
cubic box of size L and volume V = L d . We consider periodic boundary conditions. Let p 
denote the particle density, (3 the inverse temperature, and /i the chemical potential. The 
canonical partition function in the Feynman-Kac representation is given by 



N . 
Y W = E J\ E d Xl ...dx k dWgl ( Wl ) . . . dWZ*i h (u k ) 

k=l ' m,...,n k ?! i jDk 



niH \-n k =N 



n_! e -/SW(wj) 
n,- 



•i=i 



1 <C i<j sC k 



e r~-v-»-j; . (2.1) 



This expression is illustrated in Fig. ^ In words, we sum over the number k of closed 
trajectories and over their respective winding numbers m,... , n*;. We integrate over the 
initial positions x\, . . . , We integrate over trajectories Uj : [0, Uj(3] — ► D that start and 
end at xj; here, W£ x denotes the Wiener measure. See Appendix[X]for more information, 
and in particular Eq. (|A.10|) for the normalization condition. Trajectories wind around 
the time direction according to their winding numbers; because of periodic boundary 
conditions, they may also wind around space directions. 

Given a trajectory uj with winding number n, the function U(uj) denotes the interactions 
between different legs; explicitly, 



U ^)= E \l U{u>{if} + s) 
<i<i <n-l P Jo 



Lu(j(3 + s))d 



s. 



(2.2) 
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And U(uj,uj') denotes the interactions between closed trajectories u> and a/, of respective 
winding numbers n and n': 

_ i rP 

U(u,u')= V -/ U(cj(il3 + s)-u , (jf3 + s))ds. (2.3) 

The function f/(x) represents the pair interaction potential between two particles separated 
by a distance \x\. We suppose that U{x) is nonnegative and spherically symmetric. We 
can allow the value +oo; all that is needed is that e _ ^ w ^ and q - ^^'^') be measurable 
functions with respect to the Wiener measure — any piecewise continuous function D — > 
[0, oo] can be considered at this point. 
The grand-canonical partition function is 

Z(p) = ^ ^ N Y{N) (2.4) 

N ^ 

(with the understanding that Y(0) = 1). We also need partition functions where a given 
trajectory ujq is present — these will be needed in the expression for cycle densities, see 
(|2~H|) and (JZSJ). Namely, we define 



N 

Y(N; = £ d Xl ...dx k W%& K) . . . dW?*£ (w fc ) 

fc=l ' ni,...,n fc > 1 ^ 



niH \-n k =N 



k 



^ i<j ^ 



The dependence on ujq comes from the last term, where the product includes terms with 
i = 0. Notice that 

Y(N-co) ^Y(N), (2.6) 
with equality iff U(x) = 0, i.e. in absence of interactions. Finally, we introduce 

Z(fx,u)= ^ N Y{N,uj) (2.7) 

N >0 

(we set y(0,u;o) = 1). We also have Z(p,uj) ^ Z(p), with equality iff U(x) = 0. 



2.2. Cycle lengths. We now introduce the density of particles in cycles of length n, both 
in the canonical and grand-canonical ensembles. We denote the particle density by p = y. 
When discussing the canonical ensemble, we always suppose that p and V are such that 
N = pV is an integer. The number of particles in cycles of length n is given by the random 
variable Ylj=i n< ^n,-,n- Averaging over all configurations of space-time closed trajectories, 
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we get 
g p (n) 



: Y (N)^h\ E T/ E n ^ 



k\ 

k=l ni,...,rtfe ^ 1 
niH \-n k =N 



3=1 



dx\ . . . dxk 



D k 



k 

dw^(^)...dw^> fc ) 



■3=1 3 



I 



n •- 

1 < i<j < k 



<H»- w - i r( ;i„7' 



(2i 



The last line follows from the first line by replacing ^j=i n $nj,n with nkS ni:n ; isolating the 
integral over wi; using the definition (|2,5|) : using translation invariance, and J D dxi = V. 
Similarly, we have the grand-canonical expression 

-pU{w) z (p^) 



Qn{n) 

One easily checks that 



/*" d<(,)e 



z( M ) ■ 



(2.9) 



E ^ 

t ^ l 

E Sm(») = (v> = p(m). 



n > 1 



We consider the thermodynamic limits of g p (n) and Q^n). Since ^ Q P {n) ^ p, and 
since n is a discrete index, Cantor diagonal process yields the existence of a sequence of 
increasing volumes Vk, with pV^ = an integer, such that g p (n) converges to some limit 
that we denote Q p (n). Similarly, we also obtain the infinite volume limit Q^n). Fatou's 
lemma implies that 

E e » ^ p> E ^» ^ ( 2 - 10 ) 

n 1 n 1 

This suggests to define the density of particles in infinite cycles by 



Q P {°°) = P~ E q p 

n > 1 

e^(°o) = p(m) - E ^(")- 



(2.11) 



n > 1 



The main question is whether g(oo) differs from zero at given temperature, and at given 
density or chemical potential. 

We chose to discuss densities of particles in cycles of given length, but one may consider 
probabilities as well. Namely, we could introduce the probability for particle 1 to belong 
to a cycle of length n; it is given by 

NY(N) ' y ' 

Thus Q p {n) = pP p (n) in the canonical ensemble, and Q p (oo) = /xP p (oo). Things are not 
so simple in the grand-canonical ensemble. The probability P p (n) is 



P P {n) 



D 



dx / dW£f(o;)e- 



P» = dx dW£i ( W ) 



D 



^f3p,n 



n 



Z(p) 



(2.13) 
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Here, Z'(fj,;u) is like Z{\i\uS) given in Eqs (|2.7j) and (|2,5j) . but with a factor m rrn instead 
of ^. Heuristically, we should have (^r) = l/p(n), and £ p ( 



(fc+1)! 

(nj = p(fi)P^(n), but this does 
not seem easy to establish. The ratio of partition functions in (|2.13|) is more difficult to 
control than the one in (|2.9j) . We therefore abandon probabilities and discuss densities. 



2.3. Off-diagonal long-range order. Let us turn to Penrose and Onsager off-diagonal 
long-range order. Its Feynman-Kac representation involves an open trajectory that starts 
at x and ends at y, that possibly winds several times around the time direction. Precisely, 
we introduce 



o p {x) 



N 

E 

n=l 

E 

n > 1 



dW2f(a/)e 



Y(N 



n; uj 



dW2f( W )e 



y(7V) ' 



(2.14) 



Z( M ) ' 



Thermodynamic limits are denoted cr p (x) and cr^(x), provided they exist! One may 
actually restrict a p {x) and <7 p (x) on rational x, and use the Cantor diagonal process to get 
convergence on a subsequence of increasing volumes. This is not necessary in this paper, 
as the limits will be shown to exist in the regimes of parameters under consideration. 

Similarities between Eqs l|2.8|) . (|2.9|) on the one hand, and Eqs (|2.14|) on the other hand, 
are manifest. We can write 



A' 



n=l 

n> i 

where the coefficients c njP , c n ^ are given by 



n 



(2.15) 



c n,p{ x ) 



» 



| dWSf(«)e 



-/JW(u;) z (w^) 
Z( M ) 



/dO)e 



n; cj 



E(iV) 



ZGu) 



(2.16) 



As above, we denote the thermodynamic limits by c„ jP (x) and c n ^(x), provided they exist. 
One should be careful when sending the volume to infinity in Eqs (|2.15|) . because a "leak 
to infinity" may yield a term involving g(oo) — this actually occurs in the ideal gas, as 
shown in the next section. 



3. The ideal gas 

The ideal gas of quantum bosons is fascinating. Particles do not interact, yet they 
manage to display a phase transition. Historically, the Bose-Einstein condensation is the 
first theoretical description of a phase transition. The ideal gas has been the object of 
many studies over the years; let us mention |17[ QH]. A simple proof of macroscopic 
occupation of the zero Fourier mode is presented in Appendix iBl 

In this section we elucidate the relation between cycle lengths and off-diagonal long- 
range order, thus clarifying results that were previously obtained by Siito ^JE3- We work 
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in the canonical ensemble and establish the formula explicitely, for any dimension 

O l. 

Theorem 1. For any < (3, p < oo, there exists a sequence of increasing cubes for which 
the thermodynamic limits of o~ p (x), c n , p (x), g p (n) exist for all x G K d and n G N. Further, 
we have 



n > 1 



The rest of this section is devoted to the proof of Theorem ^ The coefficient c ra)P (x), 
defined in Eq. (|2.16|) . has a simpler expression in absence of interactions. Indeed, we have 
U(lo) = and Y(N;uj) = Y(N). It follows from properties of the Wiener measure in 
periodic boxes, see Eq. (|A.10|) . that 



{f-zf 



zeT, d 



(3.1) 



Notice that lim^oo CnJx) = e *™0 , but the limit is not uniform in n. If the sum over n 
is restricted to n ^ cL 2 , with c any finite constant, we can use the dominated convergence 
theorem and we get 

cL 2 _ 2 

lim Vc njP (x)e p (n) = V g (n). (3.2) 

L— >oo ' — * *■ — » 

n=l n > 1 



(The limit is taken along the subsequence of increasing volumes for which g p (n) is known 
to converge for any n.) 

We now consider the terms with cl? < n ^ N. We estimate the sums in (|3.1(l using 
integrals; we have 



/OO fOO 
-oo fegZ J-oo 



The Gaussian integral is equal to y/ir/a. Consequently, 

\y4wn(3 — L 1 d 



V4vrn/3 - L 

These bounds hold provided ^Airnf3 > L. Since > c, we have 



' yggj- 1 
We obtain 



A 



J] 6 P {n) < ^ c n , p (x)£ p (n) < 

n=cL 2 n=cL 2 



V4^b~P + 1 



v 7 !^/?- 1 



d A 



n=cL 2 



V 



cL 2 



y/4^3 + 1 

v/iTTC^- 1 



cL 2 v 
n=l ' 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



Using (|3.2() with x = and the definition Q2.11JI of the density of infinite cycles, we see 
that the last term converges to f? p (oo) as L — > oo. It then follows from (|3.2|) and (|3.6|) 
that 

lim sup o" p (£j ^ >^ e 4n P £> p (n) + 

L^oo ^ , 

n > 1 



^4^3-1 



ffp(oo). 



(3.7) 



s 



DANIEL UELTSCHI 



This inequality holds for any c, and the fraction is arbitrarily close to 1 by taking c large. 
A lower bound can be derived in a similar fashion, and we obtain the formula stated in 
Theorem ^ 



4. The interacting gas 

The interacting gas is much more difficult to study. We prove in this section the absence 
of infinite cycles when the chemical potential is negative (Theorem |2J) . We then study the 
coefficients c n ^(x) at low density and high temperature, using cluster expansion tech- 
niques. Their thermodynamic limit can be established, and we show that c n ^(x) — > as 
\x\ —* oo (Theorem |3J)- 

Theorem 2. Let < j3 < oo and \i < 0; then 
and 

lim limsupo" /J (x) = 0. 

|x|-»oo L^oo 

Proof. Since U(u) and Z(p,\u) ^ Z(/j,), the finite volume density Qp(n), Eq. ()2.9j) . is 
less than 

dw ™ {uj) = 5 e ~~ • (41) 

The right side is smaller than e^ n for all L large enough. We can therefore apply the 
dominated convergence theorem and we obtain 

p= lim V Qfl (n) = V Q (n). (4.2) 

n 1 n 1 

It follows that g p (oo) = 0. The statement about absence of off-diagonal long-range order 
can be treated similarly. We have the upper bound 

/3/m ( X -Lz) 2 

^) < E ti^w^ E e " ^ • ( 4 - 3 ) 



„Yi (47m/3)d/2 



ze 



By dominated convergence, 



lim sup ^(x) > ; 77^-^775 e^ 4 ^ . (4.4) 



e -x 2 /Anf3 



We can again use the dominated convergence theorem for the limit \x\ — > oo, and we get 
the claim. □ 

We continue the study of the interacting gas in the regime where cluster expansion con- 
verges. We assume that the chemical potential is negative, that the interaction potential 
U(x) is integrable, and that the temperature is high enough. The condition in Theorem |3] 
is stronger than necessary, but it is very explicit. We will invoke a weaker condition in the 
proof of the theorem that is based on the "Kotecky-Preiss criterion" for the convergence 
of cluster expansion. Notice that Ginibre's survey |Hj uses Kirkwood-Salzburg equations; 
it applies to a broader range of potentials, but things are terribly intricate. 
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Theorem 3. Assume that (3, \i, and U satisfy 



The thermodynamic limits of c ri ^(x) and Qfi(n) exist, and we have 



lim Cn^ix) = 

\x\— >oo 



for any n. 



Proof. We need some notation in order to cast the grand-canonical partition function in 
a form suitable for the cluster expansion. Let us introduce a measure for trajectories 
that wind arbitrarily many times around the time direction. Namely, let X n denote the 
measure space of continuous trajectories u : [0, nj3] — ► D, and let X = U„ ^ \X n be the set 
of trajectories in D with arbitrary winding numbers. We introduce the measure v on X 
whose integral means the following: 

/j3 fJ,71i p f* 
F(u)du(u) = ~ / dx / d^f(w)e-^M F(u). (4.5) 
n ^ 1 n J 

It is clear that v is a genuine measure on a reasonable measure space. But we describe 
the measure v with more details for readers who are interested in analytic technicalities. 
The cr-algebra on X n is the smallest cr-algebra that contains the sets {ui G X n : u)(t) G B}, 
for any ^ t ^ n/3, and any Borel set B C D. Trajectories of X\ can be dilated in the 
time direction so as to yield trajectories with arbitrary winding numbers. One can then 
consider the product space X\ x N with the product u-algebra (the u-algebra on N being 
the power set). The measure of a set of the kind A x {n}, with A a measurable subset of 
Xi, is defined as 

v(A x {n}) = / dx / dW"£f(a/) e-^ W(£J,) . (4.6) 



n 

Here, we introduced 

A' = {J G X n : J(t) = u(nt) for some lo € A}. (4.7) 

There is a unique extension to a measure on X\ x N. There is a natural correspondence 
between X and X\ x N, and we consider v to be a measure on X . 

With this notation, the grand-canonical partition function (|2.4j) is given by 

Z{p)= Y, h I M"i)---M"k) II [e-M^-l]. (4.8) 
k > o ' Jxh 1 < i<j < k 

The term k = is equal to 1 by definition. Then Z(\i) has exactly the form assumed e.g. 
in ^3] . The Kotecky-Preiss criterion for the convergence of the cluster expansion requires 
the existence of a function a : X — > IR + such that the following inequality holds for any 
u> G X: 



L 



3 a{uJ,) du(co') < a(cj). (4.9) 



Choosing a (a;) = —fi^n (with n the winding number of the trajectory u), it was shown in 
[Tl| that (|4.9|) is a consequence of the condition in Theorem [21 
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The main result of the cluster expansion is that the partition function Q4.8JI is given by 
the exponential of a convergent series. Namely, 



&v{ui) . . . dv{u k ) cp(uJi, ... ,Uk) 



(4.10) 



The combinatorial function tp(ux, ■ ■ ■ ,cJk) 1S equal to 1 if ft = 1, and is otherwise equal to 



G (i,j)eG 



(4.11) 



The sum is over connected graphs with k vertices, and the product is over edges of G. A 
proof for the relation (|4.1U|) that directly applies here can be found in |14j . 

Observe now that the partition function Z{\i\ uo) is given by an expression similar to 
(|4.Sjl . where each du{uj) is replaced by e~P u ^ 3 > u i> dv{ujj). Since U(u),Uj) is positive, the 
criterion (|4.9|) is satisfied with this new measure. It follows that Z(h\uj) has an expansion 
similar to (|4.1()j) . and we obtain the following expression for the ratio of partition functions, 



Z(h;uj) 



exp 



k > l 

It is not hard to check that 
k 



I dv(u; 1 )...dv(u; k )[l-lle^ u ^ 



ip(wi,...,w fc ) >. (4.12) 



i-ne 



k 

EC- 

k 
3=1 



j-1 

1=1 



(4.13) 



)• 



Then Equation (5) in |14j gives the necessary estimate for the exponent in (|4.12|) . namely 



]T / t di/(«i)...di/K)[X;(l- e-W^)]!^!,...,^)! < -Pfxn. (4.14) 



This bound is uniform in the size of the domain, which is important. It follows that, as 
L — ► oo, the ratio Z (fx; uj) / Z (fi) converges pointwise in [i and uj. The thermodynamic 
limits of c n ^{x) and £> M (n) then clearly exist. Further, c„ iM (x) is bounded by 



dW^(w)e-^H 



(4.15) 



It is not hard to show that the bracket is bounded away from zero uniformly in L (but 
not uniformly in (3 and n). From (|A.10p . we have 



lim lim / dW^{u) = 0. 



This implies that c n ^(x) vanishes in the limit of infinite |x| 



□ 
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5. Conclusion 

We introduced the formula Ijl.lJI that relates the off-diagonal correlation function and 
the densities of cycles of given length. This formula involves coefficients c n that have 
a natural definition in terms of integrals of Wiener trajectories. We conjectured several 
properties for the coefficients — these properties can actually be proved in the ideal gas 
for all temperatures, and in the interacting gas for high temperatures. These results seem 
to indicate that the order parameters of Feynman and Penrose-Onsager agree. However, 
heuristic considerations based on the present framework ^3] suggest that, if the gas is in 
a crystalline phase, the coefficients satisfy c n (x) ^ e ~ a M for some a > 0, and for all n 
(including n = oo). Besides, one expects that £?(oo) > if the temperature is sufficiently 
low. The order parameters are not equivalent in this case. 

An open problem is to establish the equivalence of the order parameters in weakly in- 
teracting gases in presence of Bose-Einstein condensation. Another question is whether 
Coo( x ) converges, as \x\ — > oo, to a number that is strictly between and 1. The corre- 
sponding phase would display a Bose condensate whose density is less than the density of 
particles in infinite cycles. 

Acknowledgments: It is a pleasure to thank Valentin Zagrebnov and Joe Pule for 
encouragements, and for pointing out an error with thermodynamic potentials in the 
original vresion of the manuscript. This research was supported in part by the grant 
DMS-0601075 from the US National Science Foundation. 

Appendix A. Feynman-Kac representation of the Bose gas 

In this appendix we recall some properties of the Wiener measure, and we review the 
derivation of the Feynman-Kac representation of the partition functions and of the off- 
diagonal long-range order parameter. A complete account can be found in the excellent 
notes of Ginibre f|]; Faris wrote a useful survey [3]. 

Let D be the d-dimensional cubic box of size L and volume V = L d . We work with 
periodic boundary conditions, meaning that D is the d-dimensionial torus T^. The state 
space is the Hilbert space T~Ld,n of square-summable complex functions on D , that are 
symmetric with respect to their arguments. Let S denote the symmetric projector on 
L 2 (D N ), i.e. 

Sip(xi, . . . ,x N ) = — ^2 ^(^(l), ■ ■ ■ ,%ir(N)), (A.l) 

where x%, . . . ,xn S D and the sum is over all permutations of N elements. The state 
space for N bosons in D is therefore TId,n = SL 2 (D N ), the projection of L 2 (D N ) onto 
symmetric functions. 

The Hamiltonian of the system is the sum H = T+V of kinetic and interaction energies. 
The kinetic energy is T = — X^jLiA?' where Aj is the Laplacian for the j-th particle. 
Interactions are given by the multiplication operator V = ^ {< j ^ N U(xi — xj). 

Recall that [5 and [i denote the inverse temperature and the chemical potential, respec- 
tively. The canonical and grand-canonical partition functions are 

Y(f3,V,N) = Tr HDiN e-P H , (A.2) 

Z(p,V,fx)= ]T e^ N Y((3,V,N). (A.3) 

N ^ 
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Under the assumption that U(x) is a stable potential and that it decays faster than \x\~ d 
as \x\ — > oo, one can establish the existence of the thermodynamic potentials (see 

f(P,p)= lim -^logF(TV), (A.4) 

p(/3,/x) = lim ilogZO*). (A.5) 

Further, / and p are related by a Legendre transform, 

/(/?,p) = sup[p At -ip(/3, j u)]. (A.6) 

A 1 

This equation is useful to find / from p in the case of the ideal gas, where p can be 
computed explicitly. 

The Feynman-Kac representation allows to express e - ^ in terms of Wiener trajectories 
(Brownian motion). We briefly review the main properties of the Wiener measure. Let 
X\ be the set of continuous paths uo : [0, 0\ — ► D. Consider a function F : X\ — > 1R of the 
kind 

FH = /(w(ti),...,w(* n )), (A.7) 

where / is a bounded measurable function on D n , and < ii < • • • < t n < f3; we extend / 
on M. d by periodicity. The integral of F with respect to the Wiener measure w£ y is given 
by 

/ F(u)dW£ '(w) = V / 5t 1 (xi-x)^ ta _ tl (a;2-a?i)---5/9-t B (!/ + ^-«n) 

/(xi, . . . ,x n )dxi . . . dx n , (A.8) 



where gt is the normalized Gaussian function with mean zero and variance 2t, 



The sum over z accounts for periodic boundary conditions. A special case of QA.8JI is when 
the function F is the constant function F(uj) = 1; we get 

/ , ^ _ (x — y + Lz)^ 

dW£ y {u) = (4^)~ d / 2 £ e — . (A.10) 

Only the term z = remains in the limit L — > oo. It can be proved that such a measure 
exists and is unique jH]. The Wiener measure W^f is concentrated on Holder continuous 
trajectories (with any Holder constant less than ^) that start at x and end at y. Integration 
with respect to w£ y and Wq Q are related as follows. Define w'(t) = u>(t) — t^^-; then 

F(u)dWP(uj) = e"T I F{J)dWg {u). {All) 



The Feynman-Kac formula states that e~@ H is given by an integral operator 
We are actually dealing with bosonic particles, and it is more convenient to consider the 
operator e~^ H S that also projects onto symmetric functions. We have 

e _/3H Si(j(xi, . . . , x N ) = / K(x 1 ,...,x N ;yi,...,y N )ip(yi,...,y N )dy 1 ...dy N , (A.12) 

Jd n 
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where the kernel K is given by 



K(x 1 ,...,x N ;yi,...,y N ) = ^ J dW^ y<x) (wj) . . . dW^ Ny<N) (u N ) 

exp|-^y C/"(wi(s)-a;j-(a))dflJ. (A. 13) 
The canonical partition function is then given by 
y ( iV )= E 4i [ d Xl ...dx N [ dWg lXv(i) (u l )...dWg NXn(N) (u N ) 

ex p|"E / ^(wi(s)-«i(s))da}. (A.14) 
1 i<j Jo J 
We now group the cycles into closed trajectories, that may wind several times around 
the time direction. The number of permutations of A elements with k cycles of lengths 
ni, . . . , n k (with JA n>j = A) is 

A! 



Further, we have 

J dx 2 ... dx n J dW&>i) . . . dW^ y (Lo n )F{Lo) = J dW$(u;)F(w). (A.15) 

The trajectory ui : [0, n/3] — > D in the right side is the concatenation of u±, . . . ,u} n . The 
partition function (|A.14j) can then be rewritten in the form l|'2.1|) . 

Let us turn to Penrose and Onsager off-diagonal long-range order [S]. Given a single 
particle function ip £ L 2 (D), we define the operator N v that represents the number of 
particles in the state p. The action of this operator is given by 

N r 

(N tp ip){xi,...,x N ) / ip(x)ip(x 1: ..., ^x^ ,...,x N )<p(xj)dx. (A.16) 

J=l ^ J-th place 

It is clear that ^ A^ ^ A, and that [N V ,S] = 0. Let <po(x) = ^= denote the single 
particle ground state in absence of interactions. It is also the Fourier function with mode 
k = 0. The average occupation of the zero mode is given by 

1 



= lim — ; — -Tr-Wn „ 
*P v^ooY(N) HD ' N l V 



N Vo -/3H 



(A.17) 



We set A = pV, and the limit exists at least along a subsequence of increasing volumes. 
A criterion for Bose-Einstein condensation is that differs from zero. We can derive a 
Feynman-Kac expression for this order parameter. From (|A.12|) . ()A.13|) . and (|A.16|) . we 

have 



TrH D A e m = ( AT 1 iM / dxp(x) f dyp{y) [ dx 2 . . . dx N V 
(A - 1)! J D J D Jd n - 1 Tt 

/ d^^O-.-d^^ (A.18) 
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Here, we set x\ = x, x\ = y, and xj = xj for 2 ^ j ^ N. Then can be written as 



<?p° )= J im T72 / a p( x ~ v) dxd v ( A - 19 ) 

F V^oo V z J D 2 

where 

ap{x ~ y) = nphr) (F^i)i d * 2 • • • dXN E 

7Tfci)iV 

/ d <^ (1) (^)--- dW ^ (JV) (^)exp{-^ / 17^)- «,■(.))}. (A.20) 

This expression involves an open cycle from x to y, winding n times around the time 
direction, with n = 1, . . . , N . Using the concatenation property (|A.15|) . and thanks to the 

combinatorial factor / jyZ^l , we obtain the expression (|2.14f) for o p {x — y). The system 
displays off-diagonal long-range order if cr p (x) is strictly positive, uniformly in V, x. 

Appendix B. A simple proof of macroscopic occupation in the ideal gas 

In this section, we give a proof of the macroscopic occupation of the zero mode at 
low temperature. This is usually established in the grand-canonical ensemble, using a 
chemical potential that varies with the volume and tends to zero in the thermodynamic 
limit. This approach is rather unnatural, and requires large deviation techniques to control 
the fluctuations of the number of particles. The present proof is simpler and stays within 
the canonical ensemble. 

The computation of the pressure and of the density in the grand-canonical ensemble 
can be found in any textbook dealing with quantum statistical mechanics. The chemical 
potential must be strictly negative. The infinite volume pressure is 

PCS. ri = ~PTs3 [ M 1 ' e-W-A ) dk, (B.l) 

and the density is 

p(A/i ) = _L f = 1 y ^ n - d ' 2 . (B.2) 

PVM ' W (27r) d J Rd e^ 2 "^) - 1 (47^/2^ V ; 

The limit of p(/3, fx) as fj> f is finite for d 3, and gives the critical density of the 
ideal Bose gas, p c . The graph of p((3,fj,) in three dimensions is plotted in Fig. |2] (a). Its 
Legendre transform (|A.6|) gives f((3,p), see Fig. [21(b); it is nonanalytic at p c . The value 
of a(/3) is given by 

a(/3) = lim ±p(/3, p) = lim -/(/?, p). (B.3) 

The very nature of the Bose-Einstein condensation is that the occupation number for 

k = becomes macroscopic. The average occupation of the zero mode gf \ see Eq. (jA.17|) . 
can be rewritten as 

Kp v^oo Y(N) V y ' 

V ' (n k ):N 

Here, N = Vp, and the sum is over all occupation numbers ^ 0, with indices k £ 
(jj-"L) d , such that ^fc n fc = N. The heart of Bose-Einstein condensation is the following 
result. 
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/3a(/3) 



-a(/J) 



(a) (b) 

Figure 2. The pressure and the free energy of the ideal gas in three dimensions. 

Theorem 4. For d ^ 3, the single-particle ground state is macros comically occupied if 
p > p c . More precisely, 

= max(0,p-p c ). 

Proof. It is clear that ^ 0. We now establish that ^ p — Pc- Let us introduce the 
average occupation of the mode k, 



{nk) = YlN) £ nfc6 
Thanks to the sum rule N = X^( n fc)> we have 



-/J£ fc , n k ,k> 



lim > 



(B.5) 



We can view as a random variable taking positive integer values; its expectation is 
therefore given by 

(n k ) = Prob(n fc > i), (B.6) 



where we defined 



Prob(nfc ^ i) 



i > 1 
1 

Y(N) 



E 



(B.7) 



{n k ,):N,i 

The sum is restricted to (nfc/) such that ^ n fc' = an d ^ i. The change of variable 
^ rik — i leads to 

Prob(n fc > <) = e"^ fe2 (B.8) 

The ratio of partition functions is also equal to the probability Prob(no ^ i), which is 
smaller than 1. Equations (|B.6|) and HB.8|) give a bound for the occupation numbers of all 
modes k ^ 0, namely, 



(nfc) < 



1 



(B.9) 



e^ 2 -1" 

Notice that (n^) ^ ^ ^ or ^ 7^ 0- This shows that only the zero mode can be 
macroscopically occupied (for d ^ 3) . Inserting this bound into (|B.5|) , we obtain 



p m >p _J_ lim vf— V - 



(B.10) 
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The limit converges to the expression (|B.2|) with p = 0, which is equal to p c . 

There remains to show that ^ max(0, p — p c ). From (|B.8|) with k = 0, and using 
the equivalence of ensembles, we have for any fixed a, 

lim -±-logProb(n >Va) = f((3,p)-f(P,p-a). (B.ll) 

V^co pv 

The right side of ()B.11|) is strictly negative when a > max(0,p — p c ). There exists 5 > 
such that for large enough volumes, 

Prob(n > Va) < e~ vs . (B.12) 



Let us assume that p — p c > 0; the case p — p c ^ can be treated similarly. Using (|B.6I) 
with k = 0, together with (|B.12|) . we get 



- Pr ob(n >i) + 77 ^ Prob(n ^ 



(no) 

It follows that Q^ 1 is less than any number o > p — p c , hence £>p°^ ^ p — p c . □ 
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